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1 Introduction 



We present in this note an update of the predictions for the top quark production cross- 
section at the Tevatron. These predictions are based on two complementary ingredients: 

1. the evaluation of the parton- level cross-sections, carried out in perturbative QCD with 
the inclusion of the full next-to-leading-order (NLO) matrix elements pQ, possibly 
improved with the resummation to all orders of perturbation theory of classes of large 
soft logarithms j2] Oj 

2. the proton parton densities (PDFs), which are typically extracted comparing existing 
data with NLO calculations available for the relevant processes, and extrapolated to 
the relevant region of Q 2 using the NLO evolution equations (more recently, accurate 
estimates of the exact NNLO results have also become available jlj, based on partial 
evaluations of the three-loop splitting functions). 

The numbers we present here are based on the theoretical framework introduced in [S| 
and [Oj, where the complete NLO calculation of the tt cross-section was improved with the 
resummation of leading [3j and next-to-leading |Hj soft logarithms appearing at all orders of 
perturbation theory. The introduction of resummation turns out to have only a mild impact 
on the overall rates (the effects at NLL are typically of the order of a few percent), but 
improves the stability of the predictions with respect to changes of the renormalization scales. 
While no progress has occurred since 1998 in the calculation itself, significant development 
has taken place in the determination of the PDFs. In addition to much improved data from 
HERA, from fixed-target DIS experiments at FNAL, and to the implementation of Tevatron 
jet and W production data in the fits, progress has occurred in the assessment of the true 
uncertainties associated with the global fits to these data. This work, which recently received 
considerable attention (Giele, Keller and Kosower ,7], CTEQ [HIE], MRST [TDJ, Botje P], 
Alekhin [12]), has led to sets of PDF parameterizations which should provide a meaningful 
estimate of the "1-cr" uncertainty deriving from PDFs to be associated to any calculations 
of hard processes in hadronic collisions. 

The introduction of these PDF sets "with uncertainties" relaxes the much constrained 
predictions which used to be anchored to predefined functional parametrizations, and it is 
natural to anticipate that the range of predictions for a given hard cross-section will be 
increased. 

2 Outline of the uncertainty estimate 

We shortly outline here the details of our calculation, before presenting the numerical results. 
Unless explicitly denoted as ctnlo, all of our results are obtained using the NLL-improved 
formalism of ref. 
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2.1 Scale uncertainty 



The evaluation of the purely theoretical uncertainty is based on the standard exploration of 
the cross-section dependence on the renormalization and factorization (/ip) scales used 
in the perturbative calculation. In this work, we follow the standard convention of considering 
the range m iop /2 < fi < 2m top , setting /jr = //p = /i. A justification for this choice can 
be found in [Hj, where it was shown that fi ~ m top /2 corresponds to a point of minimal 
sensitivity, providing a maximum of the cross-section in the range 0.1 < fi/m top < 10. In 
the range of mass consistent with the current data, and for the two CM energy values of 
run I and run II (1.8 and 1.96 TeV, respectively), the relative scale uncertainty at NLO is of 
the order of ±10%, independent to good approximation of y/S, m top and PDF sets. In this 
region of parameters, the maximum value is obtained for \x ~ m top /2, and the minimum for 
fi = 2m top - The inclusion of NLL resummation corrections reduces the uncertainty to the 
level of approximately ±5% 6 2 . This is the effect of very small NLL corrections to the NLO 
result for small values of where the NLO rate is largest, and bigger corrections for large 



For completeness, we also considered the possibility of varying independently the value 
of renormalization and factorization scale. These were chosen in the range 0.5 < /xr/ /ip < 2, 
with 0.5 < Hn t F/ m top < 2. We verified (see later) that within this range the results obtained 
using the choice //r = //p are not altered significantly, leading only to a small increase of the 
upper estimate. 

2.2 PDF uncertainty 

In the framework of jBllIIlEIj, PDFs with uncertainties come in sets of ripDF pairs, where 
npDF is the number of parameters used in the fits. Each pair corresponds to the fit obtained 
by varying of ±1(7 the value of the fit parameter eigenvalues, after diagonalization of the 
correlation matrix. By construction, the systematic uncertainty obtained for the observable 
O is given by: 



where Oi± is the value obtained using the PDF set corresponding to the variation of the 
ith eigenvalue within its error range. The central value of the prediciton is obtained using 
a reference PDF set, typically labelled with i = 0. We explore in this work the sets in 
the CTEQ6 |Hj parameterizations (ticteq — 20, corresponding to 40 sets, plus 1 reference 
set) and in the MRST 2002 [TU] compilation [n mrst = 15, corresponding to 30 sets, plus 
1 reference set). All sets in the CTEQ compilation have a s [Mz) =0.118, while those in 

2 This number, as well as all numerical estimates presented in this document, correspond to the choice 
A = 2, where A is the parameter introduced in |2] to parameterize the uncertainty about subleading higher 
order terms. In that paper, it was found that A = 2 gives a better estimate of the higher order uncertainties. 
A = 0, for example, would reduce the scale dependence to only ±2.5%, without changing significantly the 
central value of the resummed cross-section 
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Vs 


H = mtop/2 




A* = m to P 




fi = 2m t0 p 




(GeV) 


0~NLO 


0~res 


&NLO 


Ores 


0~NLO 




1800 


5.17 


5.19 


4.87 


5.06 


4.32 


4.69 


1960 


6.69 


6.71 


6.31 


6.56 


5.61 


6.11 



Table 1: Cross-section predictions (in pb) for the 1998 MRSR2 PDF and mtop = 
175 GeV. 

the MRST one have a s (Mz) =0.119. The CTEQ sets are labeled as follows: 6M for the 
default set, and 101-140 for the 20 ±lcr variations. The MRST sets are labeled as for the 
reference set, and 1-30 for the 15 ±lcx variations. In both cases, CTEQ and MRST, we use 
the default values of tolerances chosen by the two groups to best represent the uncertainty. 
In particular, CTEQ selects A^ 2 = 100, while MRST selects A% 2 = 50. 

In addition, we shall also consider three sets obtained by the MRST group in 2001 [13 , 
where the values of a s was frozen to ±ler from the central world average. We shall label 
these sets as A01L for the low-a s (a s (M z ) =0.117) fit [T3], A01H for the high-a s (a s (M z ) 
=0.121) fit [13J, J01 for a fit based on Tevatron jet data (a s (M z ) =0.121)[T3]. 

3 Results 

Table ^ summarizes the results obtained with the PDF sets used in 1998, when the work in 
ref. jH] appeared. The numbers agree with what appears in Table 1 of that document. 

Table 121 gives the central value and error for the CTEQ sets, for three values of the 
top mass (170, 175 and 180 GeV) and the two CM energies of interest (yS = 1800 and 
1960 GeV). We list the results obtained at the three reference values of the mass scale 
Tfj, = fi/rritop = 0.5, 1, 2. Table El provides the same information for the MRST sets. 

Figure ^ shows the contour plots of the NLL cross-section when //r and /xp are varied 
independently. The region defined by the oblique solid lines corresponds to 0.5 < (jlr/ 1 \if < 2. 
It shows that within this domain the range of NLL rates is compatible with the range obtained 
using n R = n F . 

In principle one should combine the uncertainty due to PDFs and that due to the scale 
choice in quadrature. We prefer to add them linearly, since the scale uncertainty is not really 
a systematic error in the strict sense. We therefore quote our range for the top cross-section 
as 

a(r^ = 2) - Aa PDF ( rfM = 2) < a < a(r M = 1/2) + Aa PDF ( rfl = 1/2) (2) 

The corresponding values are given in Table 0] The similar results for the MRST compilation 
are provided in Tabled 

Three comments are in order: 
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Vs 


m t0 p 




CT re/ (6M) 


Act 


1800 


170 


0.5 


6.22 


0.42 


1800 


170 


1 


6.10 


0.40 


1800 


170 
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5.66 


0.37 


1800 


175 


0.5 


5.29 


0.35 


1800 


175 
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5.19 


0.33 


1800 


175 
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4.81 


0.31 


1800 


180 
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4.52 


0.29 


1800 


180 


1 


4.43 


0.28 


1800 


180 


2 


4.11 


0.26 


1960 


170 


0.5 


7.97 


0.57 


1960 


170 


1 


7.83 


0.54 


1960 


170 


2 


7.29 


0.49 


1960 


175 


0.5 


6.82 


0.47 


1960 


175 


1 


6.70 


0.45 


1960 


175 


2 


6.23 


0.42 


1960 


180 


0.5 


5.86 


0.40 


1960 


180 


1 


5.75 


0.38 


1960 


180 


2 


5.35 


0.35 



Table 2: Range of cross-section predictions (in pb) for the CTEQ6 family of PDFs 
at a fixed scale r M = n/m t0 p- cr re f refers to the central value, using the 6M set, and 
Act is the error, as defined in eq. 
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m top 




°Ve/(0) 
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1960 
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1960 
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5.89 


0.19 


1960 
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5.79 


0.18 


1960 
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2 


5.38 


0.16 



Table 3: Range of cross-section predictions (in pb) for the MRST family of PDFs 
at a fixed scale r M = Li/mtop- &ref refers to the central value, using the set, and 
Aa is the error, as defined in eq. 





mtop 




<7re/(6M) 


@max 


1800 


170 


5.29 


6.10 


6.63 


1800 


175 


4.51 


5.19 


5.64 


1800 


180 


3.85 


4.43 


4.81 


1960 


170 


6.79 


7.83 


8.54 


1960 


175 


5.82 


6.70 


7.30 


1960 


180 


5.00 


5.75 


6.25 



Table 4: Full range of cross-section predictions (in pb) for the CTEQ6 family of 
PDFs, as defined in eq. ©• cr re f refers to the choice of 6M and li = mtop- 
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Figure 1: Contour plot of the NLL cross-section, in the (jlf — fJ-R plane. The oblique 
solid line defines the region 0.5 < [ar/ hf < 2. 



Vs 


mop 


@min 


cr re/ (0) 




1800 


170 


5.52 


6.13 


6.44 


1800 


175 


4.69 


5.21 


5.47 


1800 


180 


4.00 


4.44 


4.67 


1960 


170 


7.11 


7.90 


8.31 


1960 


175 


6.08 


6.76 


7.10 


1960 


180 


5.21 


5.79 


6.08 



Table 5: Full range of cross-section predictions (in pb) for the MRST family of 
PDFs, as defined in eq. ©. a re f refers to the choice of set and /i = m top . 
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(T™,v, (t,, = 2 AOIL) 


<t^ c /( r, ( = 1. 0) 

^ TG.J \ ' (1 - 1 - 5 ^ J 


Gmnrr- (T =05 JOl) 

^ max y fi km \j -l j 


1800 


170 


5.48 


6.13 


6.72 


1800 


175 


4.66 


5.21 


5.71 


1800 


180 


3.98 


4.44 


4.86 


1960 


170 


7.04 


7.90 


8.69 


1960 


175 


6.03 


6.76 


7.41 


1960 


180 


5.17 


5.79 


6.34 



Table 6: Full range of cross-section predictions (in pb) for the MRST family of 
PDFs. a re f refers to the choice of and fj, = mt op - f> = iJ./m t0 p and PDF give the 
scale factor and PDF set at which the minimum and maximum rates are attained. 



1. the uncertainty ranges obtained using the CTEQ sets, for a fixed choice of scale, are 
almost twice as large as those for the MRST sets. We understand this is the result 
of the different tolerance criteria used by the two groups (see Appendix B4 of [5] and 
Section 6 of jTU] for some discussion). The MRST range increases however if we include 
in the analsys the 2001 sets with varying a s . This is shown in Table El In this case 
the lowest predictions are obtained from the 2001 A01L fit, with the low value of a s 
, while the highest prediction comes from the 2001 jet-based J01 fit. After the a s 
variation is included, the MRST range becomes compatible with that of CTEQ's. 

2. the central values obtained today for the top cross-section are about 3% larger than 
those obtained in 1998. At y/S = 1.8 TeV and /i = m top = 175 GeV we had 5.06 pb 
with the set MRSR2 (a s (M z ) =0.119). We now have 5.19 pb with CTEQ6M, and 
5.21 pb with MRSTO. 

3. the contribution of the PDF systematics to the uncertainty range is large. In the case 
of the CTEQ sets, it is of the order of 6-7%, larger than that due to the choice of scale. 
This is a result of the large sensitivity of the top cross-section to the large-x gluon 
content of the proton, which is still poorly known. For CTEQ the largest contribution 
to the error comes from the two sets 129 and 130 3 . For these two sets, we find the 
contribution of the gg channel to be respectively 11% and 21% of the total rate. For 
comparison, the contributions of the qq production channel for sets 129 and 130 are 
the same to within 1%. In other words, the PDF uncertainty on the top rate is mostly 
driven by the poorly known gluon density, whose luminosity in this kinematic range 
varies by up to a factor of 2 within the 1-a PDF range. 

While the overall production rate has a large relative uncertainty of approximately ±15%, it 
is important to point out that the ratio of cross-sections at y/S = 1.96 TeV and y/S = 1.8 TeV 
is extremely stable. In the case of the CTEQ sets, for example, we found <r(1.96)/cr(1.8) = 
1.295 ±0.015 after scanning over the set of scale choices and for 170 < m top < 180 GeV. The 

3 This is consistent with what found in a recent study of jet produciton at the Tevatron |14j . 
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' ' v top 




a rP f (6M.) 


^vnax 


1800 


170 


5.29 


6.10 


6.72 


1800 


175 


4.51 


5.19 


5.71 


1800 


180 


3.85 


4.43 


4.86 


1960 


170 


6.79 


7.83 


8.69 


1960 


175 


5.82 


6.70 


7.41 


1960 


180 


5.00 


5.75 


6.34 



Table 7: Full range of cross-section predictions (in pb) for the combined study of 
CTEQ6, MRST and MRST with a s variation. The central vlaues are taken from 
CTEQ6M. The minimum rates arise from CTEQ6, while the upper values arise from 
MRST set J01. These numbers should be quoted as "BCMN as updated in [this 
paper]." 

error is about 1%. We therefore consider the prediction of the relative cross- section at the 
two energies to be a very stable one. 

For reference, we collect the full set of cross-sections (at y/S = 1.96 TeV and m top = 
175 GeV) for all CTEQ sets and scale choices in Table|Hl Here, for the sake of documentation, 
we provide the NLO rates and the NLL-improved ones separately. 

4 Conclusions 

We reiterate here the main findings of this study. The inclusion of the full PDF systematics, 
made possible by the recent works of several groups, leads to a more realistic estimate of the 
top cross-section uncertainty. The latest MRST and CTEQ sets give rise to cross-sections 
which are typically 3% larger than what obtained with sets available at the time of Run I. 
In addition to the increase in rate, the size of the uncertainty range has also increased, to 
a value of the order of ±15%, dominated by the PDF and a s uncertainties. The leading 
source of PDF uncertainty comes from the (lack of) knowledge of the gluon luminosity at 
large values of x. The gg contribution can in fact change through the PDF sets by up to a 
factor of 2 (from 10% to 20% of the total rate at 1.96 TeV). We find that the MRST sets give 
rise to a smaller PDF uncertainty, a result we ascribe to the tighter tolerances required by 
MRST in defining the range of the eigenvalues. The MRST uncertainty increases however 
to values consistent with CTEQ's once the sets obtained from a ±lcx change of a s (Mz) are 
included. This underscores the importance of including the a s uncertainty into the PDF fits 
in a more systematic fashion. On the same footing, the impact of higher order corrections, 
as well as of the treatment of higher twist effects in the fitting of low-Q 2 data, may need 
some more study before a final tabulation of the PDF uncertainties is achieved • 

We collect in Table [7| our final results. This summary table includes the CTEQ6M set 
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and /i = rritop as central values, and the most extreme rates extracted from Tables HI El and 
El as lower (with fi = 2m top ) and upper values (with \i = m top /2). 

In spite of the overall large uncertainty, the ratio of cross-sections at 1.96 and 1.8 TeV is 
extremely stable, being equal to 1.295 ± 0.015 over the mass range 170 < m top < 180 GeV. 
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6.68 


5.75 


6.21 



Table 8: Full set of predictions for the CTEQ family of PDFs, and for mt op = 
175 GeV, at = 1.96 TeV. ctnlo is the NLO rate, while a res is the sum of NLO 
and NLL resummed, according to 6 . All rates are in pb. 
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